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Motivated by a recent neutron scattering experiment on detwinned YBa2Cu306+a superconduc- 
tor, we examine the frequency and doping dependence of the anisotropy in the spin and charge 
fluctuation arising from the coupling between the plane and the chain. Starting from the two- 
dimensional t-t - J model and using the random-phase approximation (RPA), we find a pronounced 
anisotropy of the incommensurate (IC) peaks in the spin channel, namely the peak intensity at the 
(tt ± 8, 7r) direction is stronger than that at the (n, n ± 8) direction in a wide frequency range from 
uj = 0.2J to the resonance frequency ui r = 0.35J. Above the resonance frequency, the IC peaks 
reemerge. Their intensities shift to the diagonal direction and no anisotropy exists along the two 
diagonal directions. We find that this anisotropy is robust with respect to the possible variation of 
the RPA correction factor and to the dopings. The charge fluctuation is also found to be incommen- 
surate for all energies considered and peak at (0, 8) and (8, 0). An anisotropy in its IC peak intensity 
along the k x and the k y direction exists, but in sharp contrast to the spin channel, the maximum 
intensity of the IC peak is along the k y direction. Moreover, the IC peak in the charge channel 
exhibits an upward dispersion, in contrast to the downward dispersion below the spin resonance 
frequency for the spin IC peak. We explain these features based on the effect of the plane-chain 
coupling on the topology of the Fermi surface. 

PACS numbers: 74.72.Bk, 74.25. Ha, 74.20.Mn 



I. INTRODUCTION 

The spin and charge dynamics are important issues 
in the studies of the high-T c superconductivity. The in- 
elastic neutron scattering (INS) experiments have been 
playing important roles on studying their spin dynamics. 
It can provide direct information of the momentum and 
frequency dependence of the dynamical spin susceptibil- 
ity. On the other hand, the INS experiments can measure 
the charge fluctuations indirectly by observing a change 
in the mass density. In addition to the indirect measure- 
ment, the charge susceptibility can be detected directly 
by the inelastic x-ray or electron scattering. 

Over the past decade, the INS experiments carried out 
on the twinned samples have observed the four-fold sym- 
metric incommensurate (IC) peaks at the wave vector 
q = (tt, 7t ± 5) and q = (tt ± 5, tt) in a range of low 
energies 0, 0, , m addition to a sharp commensurate 
resonance peak around the antiferromagnetic (AF) wave 
vector Q = (tt, tt) at about 41 meV 0,|f||(J- The experi- 
ment also reported that the IC peaks at (0, ±25) are ob- 
served for the charge excitation through the measurement 
of atomic mass density change 0,H|. Theoretically, there 
are two possible explanations of the spin incommensura- 
tion in the INS experiments. One is based on the striped- 
phase picture, which suggests that the presence of the 
dynamic stripes is the origin of the IC peaks 0, . The 
other is the scenario of the nested-Fermi-surface which 
can explain the gross feature of the incommensurate spin 
fluctuations and spin resonance [TllO ll l ll l ll MIMIr7t . 



Recently, Mook et al. have carried out the neu- 
tron scattering measurement on partly detwinned 
YBa2Cu306.6 UM an d found that the IC peak shows the 
one dimensional (ID) feature. They claimed that it gives 
a strong support on the scenario of the stripe-phase pic- 
ture. However, we showed that [l^ an anisotropic IC 
peak can in fact be obtained when the coupling between 
the plane and chain, and the asymmetry of the crystal 
lattice between the a and b axis are considered in the 
two-dimensional t-t - J model, based on the Fermi surface 
topology. Very recently, an INS experiment on fully de- 
twinned YBCO 20] shows that the geometry of the mag- 
netic fluctuations is actually two-dimensional, but with a 
prominent anisotropy in the amplitude and the width of 
the IC peak which is maximum along the (q, tt) direction 
and the anisotropy is strongly energy-dependent. Our 
previous theoretical calculation [l^ shows many features 
which are consistent with the experimental data |20| . 
However, the linewidth anisotropy found in the exper- 
iment is more stronger than our calculation at the fixed 
energy oj — 0.3J(J is the AF coupling constant in the 
t-t - J model) and a detail theoretical investigation of the 
energy dependence of the linewidth anisotropy as found 
experimentally still lacks [21| . Another new development 
in recent INS experiments is the reemergence of the IC 
peak at higher energies above the spin resonance en- 
ergy H3 ■ This high-energy IC peak is found to shift 
from the (tt, g)(and its symmetric points) direction to the 
diagonal direction |22j, |23| . This reemergence of the spin 
incommensuration has in fact been predicted theoreti- 
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cally based on the scenario of the Fermi surface topol- 
ogy [l3l IT3. fl5| . However, both the experimental and 
theoretical investigations are carried out in the twinned 
systems. So, an interesting issue is how the IC peak 
behaves in this high-energy region in the detwinned sys- 
tems. 

Though we have shown before that the plane-chain 
coupling, and the asymmetry of the crystal structure 
can lead to the in-plane anisotro py o f the spin IC peak, 
their effect may differ in detail [lJj. Namely, the for- 
mer leads the maximum IC peak to be along the k x di- 
rection, but the latter be along the k y direction. Due 
to the much stronger anisotropy caused by the plane- 
chain coupling, the combined effect gives rise to a maxi- 
mum IC peak alon g th e k x direction, which is consistent 
with experiments \18L l20j . Our investigation is based 
on a self-consistent calculation of the quasiparticle dis- 
persion and the superconducting (SC) gap in the frame- 
work of the slave-boson approach to the t-t - J model. We 
note that a recent calculation based on the similar RPA- 
type theory with the phenomenological quasiparticle dis- 
persion and SC gap obtained by fitting to experiments 
shows that the asymmetric crystal structure also leads 
to the fc^-direction IC peak [24] . in contrast to what we 
obtained before. We find that this contrasting behav- 
ior comes from the different relative ratios between the 
a — b anisotropy of the nearest-neighbor hopping constant 
and the SC gap asymmetry, both are caused by the or- 
thorhombic crystal structure. The asymmetric SC gap 
leads to a fcj,-direction IC peak, while the anisotropy in 
the hopping constant leads to a fc^-direction IC peak. 
Our self-consistent calculation gives a A x / A y w 1.5 (the 
gap function is Ak = A^ cos k x — A y cos ky), which is 
bigger than the value used by Eremin et al. [2J]. So, the 
overwhelming asymmetric SC gap in our calculation leads 
to the anisotropy be in the k y direction. Therefore, the 
explanation of the anisotropy in the IC peak based on 
the orthorhombic crystal structure may depend on the 
fine tuning of the dispersion and SC gap, in particular 
its effect is shown to be much weaker than that from the 
plane-chain coupling [l^. So, we will just focus on the 
effect of the plane-chain coupling in this paper. 

In this paper, we will first extend our previous 
scheme to address the above issues in view of the new 
experimental results pol I22L |23| . We will also examine in 
detail the effects on the anisotropic IC peaks of dopings 
and the RPA correction factor a(see the text below) [11] . 
in order to check the robustness of this anisotropy. Our 
aim is not only to account for the new findings observed 
by the very recent experiment |20|, but also to predict 
some new features which help to determine if the scenario 
of the nested-Fermi-surface can describe all features of 
the spin incommensuration or some exotic scenarios such 
as the stripe phase or the electronic nematic phase [2^ 
may be required. In addition, we will also calculate the 
charge susceptibility and investigate the frequency de- 
pendence of the anisotropy in the charge channel. 

The article is organized as follows. In Sec. II, we intro- 



duce the model and define notations. The Hamiltonian 
and formalism are obtained. In Sec. Ill, we present nu- 
merical results of the spin and charge susceptibility, re- 
spectively. In Sec. IV, we interpret the results. Finally, 
we present the conclusion in Sec. V. 

II. HAMILTONIAN AND FORMALISM 

We start with a Hamiltonian which describes the sys - 
tem with a Cu02 plane and a CuO chain per unit cell ^3 • 

H = H t _ t ,_ J +H c + H I , (1) 

where H t _ t > _j is the two-dimensional t-t - J Hamiltonian 
describing the Cu02 plane, 

Ht-t'-J = ~t'^ / c ia C i+ao- — h.C. — t c\ a Cj+ x + ya ~ h.C 

iota ia 

+J J^(Si • S l+a - -niiii +a ), (2) 

ia 

H c describes the CuO chain, 

H c = t c d\ a d i+a(X - h.c, (3) 

icta 

and the coupling between the plane and chain is p6| , 
Hi = -t±J2(cld ja + h.c) - A/4^(A f M Q A 2/J+E/2iQ 

+A t llQ A 2il _ e/2:Ct + h.c), (4) 

where c\ a and d\ a are the creation operators of elec- 
trons in the plane and chain, respectively, a — x,y 
stands for unit lattice vector along the a and b axis, 
and Si = \c\ a (j a fiCip . A nra are the singlet pair opera- 
tors, which are expressed by Ai ra — c r |C r + Q | — c r |C r+Q | 
and A 2rQ = d r ^d r+a i — d,.|d r+Q |. Following Ref. [T^ |. 
we use the slave-boson mean-field theory to the plane 
Hamiltonian H t _ t >_j, in which the physical electron 
operators ci G are expressed by slave bosons bi carry- 
ing the charge and fermions fi a representing the spin, 
Cia = blfi a . The mean-field Hamiltonian can be writ- 
ten in a matrix form [13 . 127] , and the mean-field order 
parameters A uj =< fiffn - Ai/iT >= ±^i- A 2ij =< 
di-[dji~diidj-\ >— ±A 2 (± depend on if the bond < ij > 
is in the x or the y direction), Xij = S CT < f} a fj a >= Xo, 
and the chemical potentials fi± and /i 2 for different pla- 
nar doping densities x can be obtained from the self- 
consistent equations The doping densities for the 
chain are fixed to x c h — 0.5. The other parameters are 
choosed as t = 2J, t = — 0.45t, A = 3J, the renormal- 
ized plane-chain hoping constant t± = tj_^/x = — 0.1 J, 
and the temperature T = 0.005 J. The chain is along k y 
direction. 
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The bare spin susceptibility Xg(q, c<j) coming from 
the particle-hole excitations of the planar fermions can 
be calculated from the two-particle Green's function 
Xo(Qi by the analytical continuation iuj —> uj + iS, 

k i j 

+U U (k) U 2i (k) U-Lj (k + q) [/ 2j (k + q)] 
/(^(k + q))-/^)) 
w-E i (k + q)+E i (k)+i<J' 1 ' 

with C/y are the elements of the matrix which diagonal- 
izes the mean- field Hamiltonian |l9| . 

The correction of the AF spin fluctuations to the spin 
susceptibility is included in the form of the renormalized 
random-phase approximation (RPA), 

X s (q, w) = Xo(q> w)/[l + a J q X§(q, w)]. (6) 

Here J q = J(cosg x + cosq y ). The RPA correction factor 
a (a < 1) is introduced to overcome the overestimation 
of the AF correlation in the usual (a = 1) RPA [TTll2^ |. 
Experimentally, the spin resonance peak is observed at 
41meV at the optimal doping x — 0.14 which corresponds 
to about 0.35J. Thus, whenever not be indicated, the 
AF correction factor a — 0.5 will be chosen to make the 
resonance frequency occur at 0.35 J at the optimal doping 
density x = 0.14. 

The bare charge susceptibility can also be calculated 
from the two-particle Green's function, 

k i j 

-UhQQUvMUuQl + q)C/ 2j (k + q)] 
/(^(k + q))-/^)) 

In the form of RPA, the renormalized charge susceptibil- 
ity can be written as 

X c (q, w) = xS(q, «)/[l - \ J q xS(q, «)]. (8) 

III. THE NUMERICAL RESULTS 

A. dynamical spin susceptibility 

Fig.l(a-c) show the intensity plot of the imaginary part 
of the spin susceptibility in the momentum space with 
doping x — 0.14 at frequencies 0.3J, 0.5J, and 0.75J, re- 
spectively. The IC peak intensity as a function of fre- 
quency at different momentum directions is shown in 
Fig. 1(d). From Fig. 1(a), one can see that the IC peaks 
show clear anisotropy between the k x and k y direction 
at uj = 0.3J, namely, the incommensurability S and the 
intensity of the IC peaks at (w ± 5 X , w) is obviously larger 
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FIG. 1: Panels a)-c) are the intensity plot of the imaginary 
part of the spin susceptibility Im^ 3 as a function of the wave 
vector with doping x — 0.14 for uj — 0.3J, 0.5J, and 0.75J, 
respectively. Panel d) shows the frequency dependence of the 
maximal intensity of Imx s - The inset of panel d) shows the 
ratio of the maximal Imx s at different momentum directions 
(see text). 



than that at (tt,it ± 6 y ). When frequency is increased 
to be above the resonance frequency 0.35 J, the IC peaks 
rotate to the diagonal direction, as is shown in Fig. 1(b) 
and Fig. 1(c). Interestingly, now no anisotropy exists 
along the two diagonal directions. The anisotropy be- 
low the resonance frequency is found to be frequency 
dependent, as shown in Fig. 1(d). Detailed frequency 
dependence can be found in the inset of Fig. 1(d), in 
which the ratio of the peak amplitude between the k x 
and k y direction (771), and between the diagonal and the 
k v direction (772), as a function of energy u is presented 
(771 (w) =maxqlmx s {q,TT,uj)/ma,Xqlmx s (7T,q,uj),ri2(uj) = 
max,jlmx s (<7, q, lu) / max q lmx s (tt , q, u>)). As shown, the 
anisotropy is significant in a wide frequency range from 
slightly below the resonance frequency to about 0.2J. 
It increases firstly when the energy is reduced from the 
resonance and reaches the maximal value r\\ — \.l when 
frequency is at about 0.28J, then decreases with the fur- 
ther decrease of energy. The frequency dependence of 
the anisotropy is reasonably consistent with the exper- 
iment 20], in which the maximal ratio is about 2.1 at 
u> = 28 meV. Below 0.2J, the IC peak rotates again to 
the diagonal direction as expected from the node-to-node 
excitations at low frequencies. 

We now study the dependence of the anisotropy on the 
possible variations of the AF correction factor a with fre- 
quency uj = 0.3 J and on doping density x with frequency 
u) = 0.25J in Fig.2(a-b). The main panel shows the max- 
imal intensity of spin excitations at different momentum 
directions and the inset is the ratio of the peak amplitude 
between the q y = ir and the q x = tt direction^). From 
Fig. 2 (a), one can see that the anisotropy is enhanced 
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a) 




FIG. 2: Dependence of the maximal peak in Inix" on the RPA 
correction factor a [a)] and doping x [b)], with ui = 0.3 J. The 
inset of figures shows the ratio of the maximal Imx s between 
the (q, 7r) and the (n, q) direction. 



gradually with an increasing of the AF correction factor 
a. When a equals 0.47, the anisotropy reaches its max- 
imal value. Then, the further increase of a makes the 
anisotropy decrease and when a = 0.57, the anisotropy 
vanishes. As shown in Eq.(6), the AF correction factor a 
acts to enhance the AF correction, so the resonance peak 
is shifted to lower energies with the increase of a. When 
a = 0.57, the spin resonance right appears at u) = 0.3 J, 
so the anisotropy vanishes because the spin resonance 
peaks at (w,tt). The critical doping density for the AF 
instability observed by experiments is x c = 0.02 ~ 0.05. 
It will give rise to a = 0.34 ~ 0.43 if one determine the 
AF correction factor a according to this experimental 
value x c 0,0,0- In this range of a, the anisotropy 
is about 1.35 ~ 2.0, as can be seen from the inset of 
Fig. 2 (a), so it is robust with respect to the possible vari- 
ation of a. The doping dependence of the anisotropy 
with frequency lo — 0.25J is shown in Fig. 2(b). Re- 
markably, the anisotropy increases with doping in the 
underdoped regime, then decreases with doping and still 
have a ratio r\\ — 1.2 at the heavily overdoped regime 
x = 0.26. This specific behavior is due to the variation of 
the topology of the Fermi surface with doping. Because 
this dependence is experimentally accessible, it will be 
interesting to compare the experimental result with our 
theoretical prediction, considering that the recent exper- 
iment was just carried out at the nearly optimal doped 
YBa 2 Cu 3 6 . 85 M- 

Fig.3(a-c) show the intensity plot of the imaginary part 
of the spin susceptibility Imx s as a function of the fre- 
quency and momentum. At the resonance frequency, the 
peak position is at (tt, it). Below and above the reso- 
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FIG. 3: Intensity plot of Imx s as a function of frequency ui 
and wave vector q. Panels a)-c) are along the (n,q), (q, n) 
and (q, q) direction, respectively. 



nance frequency the IC peaks develop. The incommen- 
surability is bigger as the frequency is gradually far away 
from the resonance frequency. Therefore, they exhibits a 
downward dispersion for frequency below the resonance 
frequency and an upward dispersion for frequency above 
the resonance frequency. This feature is consistent with 
the recent neutron scattering data [22ll23| . From Fig. 3(a) 
and Fig. 3(b), we can see an obvious spin gap at the (q, it) 
and (-7T, q) direction. Compared with that calculated for 
the system without the plane-chain coupling 0] , a pro- 
nounced anisotropy in the spin gap is found, namely, the 
spin gap in the (q, n) direction is bigger than that in 
the (n, q) direction. Due to the existence of nodes, there 
is no spin gap along the diagonal direction as shown in 
Fig.3(c). 



B. dynamical charge susceptibility 

In Fig.4(a-c), we show the imaginary part of the charge 
susceptibility Im% c as a function of the momentum with 
doping x — 0.14 at frequencies 0.1J, 0.3J, and 0.5J, re- 
spectively. An obvious feature of the charge susceptibil- 
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FIG. 4: Panels a)- c) are the intensity plot of the imaginary 
part of the charge susceptibility lmx c as a function of wave 
vector with doping x — 0.14 for u) = 0.1J, 0.3 J and 0.5 J, 
respectively. Panel d) shows the frequency dependence of the 
maximal intensity of Imx c . The solid line denotes the result 
along the (0, q) direction and the dashed line along the (q, 0) 
direction. 



ity is that its IC peak appears at small momenta around 
(0, 0), namely, at (0, S) and (S, 0). However, in sharp con- 
trast to the anisotropy in the spin channel, the maximum 
IC peak in the charge channel is along the k y direction. 
Meanwhile, the IC peak is always along the (0, 6) and 
(<5,0), while that in the spin channel will rotate to the 
diagonal direction at low and high frequencies. The fre- 
quency dependence of the IC peak intensity is shown in 
Fig. 4(d). As shown, the anisotropy exists in a wide fre- 
quency range. It firstly appears at lo c — 0.1J and be- 
comes significant above u = 0.2 J. 

The intensity of the charge susceptibility as a function 
of the momentum and frequency is shown in Fig. 5(a) and 
Fig. 5(b). The charge susceptibility is incommensurate for 
all frequencies considered. The incommensurability will 
increase as the frequency increases and form a usual up- 
ward dispersion. This is quite different from the disper- 
sion found in the spin channel, where the spin resonance 
at the commensurate AF wave vector Q = (tt, 7t) disparts 
the dispersion into a downward form at low frequencies 
and an upward form at high frequencies. In the stripe- 
phase picture for the spin IC fluctuation a change in 
the AF phase of the stripe at the charge-domain bound- 
ary occurs and it leads to a spin periodicity twice that of 
the charge periodicity. As a result, a two-to-one ratio of 
the incommensurability 8 between the charge fluctuation 
and the spin fluctuation is expected. Here, though we can 
obtain the relation at a certain frequency uj = 0.3 J, gen- 
erally there should be no such correspondence because 
of the different curvature of their dispersion. This fea- 
ture comes from the topology of the Fermi surface and 
its origin will be discussed in the following section. It 
constitutes an additional difference between the Fermi- 




FIG. 5: Panels a) and b) are the intensity plot of Imx c as a 
function of frequency ui and wave vector q along (0, q) and 
(q, 0) direction, respectively. Panel c) is the frequency depen- 
dence of ImXc with different wave vectors. The inset of panel 
c) shows the frequency dependence of Imxo- 



surface explanation and the stripe-phase explanation. In 
Fig. 5(c), the imaginary part of the charge susceptibility 
as a function of the frequency is shown at small momen- 
tum transfers. For smaller momentum transfer q, the 
charge excitation exhibits a broad spectrum. With the 
increase of q, the lower energy spectrum is suppressed 
and a peak occurs at the edge of the gap. With the in- 
crease of q, the peak intensity increases and reaches the 
maximum intensity around |q| = 0.47T, as can be seen 
from Fig. 5(a) and (b). From the inset of Fig. 5(c) we can 
see that the imaginary part of the bare charge susceptibil- 
ity also has a peak at the same frequency. Moreover, we 
have checked the pole equation 1 — l/4J q Xo(q, w) = 0, 
and found that it does not satisfy at the peak energy. 
So, this peak is not the resonance peak as observed in 
the spin channel, rather it is right the IC charge peak. 
Thus, an anisotropy in the peak frequency and width ex- 
ists in two different directions of the momentum transfer 
q, as shown in Fig. 5(c). 
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FIG. 6: Contour plot of the quasiparticle energy — u/2 
shown in the first quadrant of the Brillouin zone for lo = 
0,0. 12J, 0.3J and its by (7r,7r) shifted image at ui — 0.3J. 
The dashed line shows the normal state Fermi surface. The 
solid arrows denote the nesting vectors relative to (tt, tt). The 
dotted arrows denote the nesting vectors within the energy 
contour. 
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FIG. 7: Panel a) shows the imaginary part of the bare spin 
susceptibility as a function of wave vector with ui = 0.3J. 
Panel b) is its real part. 



IV. DISCUSSION 

Let us now explain the origin of the above features 
based on the topology of the Fermi surface. In the 
Fermi-liquid based theory, the spin and charge suscep- 
tibilities arise from the scattering of the SC quasipar- 
ticle across the SC gap. Their only difference is the 
sign difference in their coherence factors, namely, c s ^ — 
1 - Ok£k+q ± A k A k+q )/£; k £; k+q , which is reflected by 
the plus/minus sign in front of Uu in Eq.(5) and (7). 
Therefore, the coherence factor is appreciable when the 
gap A k and A k+q have the opposite sign for the spin 
channel and the same sign for the charge channel. Thus, 
the peak of the imaginary part of the spin susceptibil- 
ity will be near the AF wave vector Q = (tt, tt) because 
A k = — A k +Q and that of the charge susceptibility will 
be near the (0, 0) point. Then, the origin of the observed 
IC peaks can be illustrated in term of the initial and final 
scattering processes along energy contours, with different 
transfer wave vectors for the spin and the charge channel. 
Fig. 6 shows the contour plot of the quasiparticle energy 
in the SC state E k = uj/2 at u = 0, 0.12 J, 0.3 J, and its 
by (tt, 7t) shifted image at u> = 0.3 J. In the spin channel, 
the transfer wave vector of the main scattering processes 
should be near (tt, tt). Therefore, for an intermediate ex- 
citation energy such as u = 0.3 J, the best nesting vector 
(solid arrows in Fig. 6) will be a horizontal and vertical 
offset to (tt, 7t), which connects the flat part of the energy 
contour with its by (tt, tt) shifted image. The scattering 
with this best nesting vector contributes directly to the 
IC peaks, therefore the IC peaks will appear in the par- 
allel direction. In the case that there is no plane-chain 
coupling, the energy contour will be symmetric with re- 



spect to the diagonal direction. The coupling between 
the plane and chain will break this symmetry and the 
resulting nesting portion of the Fermi surface connected 
by the nesting vector along the k x direction is larger than 
that along the k y direction, as shown in Fig. 6. So, both 
the imaginary and the real part of the bare spin sus- 
ceptibility Xo show an obvious anisotropy, i.e., the peak 
intensity at (tt±5 x , tt) is stronger than that at (tt, tt±S v ). 

However, the anisotropy coming just from Im^o is in 
fact small, as can be seen from a comparison between 
that in Fig. 7(a) and Fig. 1(d). The RPA correction fac- 
tor 1 + aJqRexo i n the renormalized spin susceptibil- 
ity Imx s acts to enhance this anisotropy. The minimum 
of the RPA correction factor along the (q, tt) direction 
is less than that along the (tt, q) direction, because the 
maximum value of Re^o along the former direction is 
larger than the latter direction, as shown in Fig7(b). As 
frequency approaches to the resonance frequency, Rex§ 
will tend to — l/(aJ q ) and this makes the IC peak move 
close to (tt, tt). Therefore, the anisotropy will become 
less and less and eventually vanish at the resonance en- 
ergy. In addition, the RPA correction factor also affects 
strongly the incommensurability S of the spin susceptibil- 
ity. From Fig. 7 (a), one can notice that the IC peak po- 
sition in Im^o along the (tt ± S Xo , tt) direction is nearer to 
the (tt, tt) point than that along the (7r, 7r± <5 yo ) direction 
(S Xo < S Vo ). While in the renormalized spin susceptibility 
[shown in Fig. 1(a)], this anisotropy of the incommensu- 
rability is reversed, i.e., S x > S y . In the experiment po| 
one can also find that the peak position at (tt, tt ± S y ) is 
a little nearer to (tt, tt) than that at (tt ± 6 X , tt). 

At low frequencies, the RPA correction factor plays 
minor role. The anisotropy of the spin gap along k x and 
k y direction is mainly determined by the anisotropy of 
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the bare spin susceptibility. As frequency tends to zero, 
only node-to-node excitation becomes available, so there 
is no spin excitation along the (q, tt) and (tt, q) direction. 
As the frequency increases to ui y — 0.12 J, the constant 
energy contour crosses the line k x — tt, as shown in Fig. 6. 
So, the spin excitation along (tt, q) appears because now 
this constant energy contour and its by (tt, tt) shift image 
[not shown in Fig. 6] can be connected by the nesting wave 
vector along the k y direction. The critical frequency u y 
may be taken approximately as the spin gap along the k y 
direction. Because the constant energy contour does not 
cross the line k y = tt at this frequency, there is still no 
spin excitation along the k x direction, i.e., the spin gap 
along the k x direction is bigger than that along the k y 
direction. 

In the high energy region, the energy contour will 
evolve into two constant energy lines which are parallel 
approximately to the Fermi surface in the normal state, 
so the best nesting wave vector will rotate to the diago- 
nal direction. In this case, the nesting wave vector along 
the (tt + 8, tt + 8) and the (tt + 8, tt — 8) direction is sym- 
metric. Thus, the plane-chain coupling can not cause the 
anisotropy in the diagonal direction in the Fermi-liquid 
based theory. 

In the charge channel, the transfer wave vector of the 
main scattering processes should be near (0,0). There- 
fore, the best nesting vector (dotted arrows in Fig. 6) is 
within the energy contour and contributes the IC peak 
at (0, 8) and (8, 0). At low energies, the plane-chain cou- 
pling is less effective, so the constant energy contour is 
basically symmetric with respect to the diagonal direc- 
tion and there is no anisotropy in the charge fluctuation. 
As the frequency increases, the plane-chain coupling de- 
structs part of the nesting portion. It leads the nesting 
portion of the Fermi surface along the k y direction to be 
larger than the k x direction, so the IC peak along the k y 
direction is stronger than that along the k x direction. 

V. CONCLUSION 

In summary, we examine the frequency evolution of 
the spin and charge susceptibility based on the two di- 



mensional t-t'-J model considering the coupling between 
the plane and chain. In the spin channel, the IC peaks 
appear at (tt, tt ± 8 y ) and (tt ±S Xl tt) below the resonance 
frequency(0.2 J < ui < 0.35 J). In this region, a pro- 
nounced anisotropy of the spin susceptibility is found, 
i.e., the peak intensity along the k x direction is stronger 
than that along the k y direction. Above the resonance 
frequency, the IC peaks reemerge and rotate to the di- 
agonal direction, and there is no anisotropy in the IC 
peak intensity along the two diagonal directions. The 
dependence of the anisotropy on the antiferromagnetic 
correction factor a and the doping density x is investi- 
gated, and it is found that the anisotropy is robust in 
the reasonable range. In addition, we also find an ob- 
vious anisotropy of the spin gap between the (tt, q) and 
the (q, tt) direction. The spin gap along the k x direction 
is bigger than that along the k y direction. In the charge 
channel, the susceptibility is also incommensurate for all 
energies we considered. The IC peaks appear at (0, ±<5) 
and (±8, 0). In sharp contrast to the spin channel, the IC 
peak intensity in the charge channel along the k y direc- 
tion is stronger than that along the k x direction. Mean- 
while, the IC peak exhibits an upward dispersion in the 
whole energy region, in contrast to the downward disper- 
sion below the resonance frequency in the spin channel. 
In addition, no resonance mode is found in the charge 
channel. We explain all the results based on the scenario 
of the nesting Fermi surface and suggest that the cou- 
pling between the plane and chain is responsible for the 
observed anisotropy. 
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